Effects of free will and massive opinion in majority rule model 
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Abstract 

We study the effects of free will and massive opinion of multi- agents in a majority rule model 
wherein the competition of the two types of opinions is taken into account. To address this 
issue, we consider two specific models (model I and model II) involving different opinion-updating 
dynamics. During the opinion-updating process, the agents either interact with their neighbors 
under a majority rule with probability 1 — g, or make their own decisions with free will (model I) 
or according to the massive opinion (model II) with probability q. We investigate the difference 
of the average numbers of the two opinions as a function of q in the steady state. We find that 
the location of the order-disorder phase transition point may be shifted according to the involved 
dynamics, giving rise to either smooth or harsh conditions to achieve an ordered state. For the 
practical case with a finite population size, we conclude that there always exists a threshold for 
q below which a full consensus phase emerges. Our analytical estimations are in good agreement 
with simulation results. 
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For a long time in history simple statistical models are used to study complex biological, 
social, and geological phenomena . Recently, a great deal of efforts are devoted to 

mathematical modelling of collective behaviors of individuals, particularly opinion spread- 
ing and formation, by using for example two-state interacting spin models 
In this context, the common formulation assumes that the agents are located at the nodes 
of a network, and are endowed with two states, i.e., spin up and spin down, which mimic 
human attitudes or decisions: like/dislike, agree/ disagree, accept /reject, etc. The interac- 
tions are assumed to take place only between linked nodes, where the links represent the 
relations between individuals such as friendship, coauthorship, neighborhood, etc. Through 
interactions, the agents can update their opinions by ways of convincement [sl, imitation Q], 



following local majority 



4, la, [3], and so on. Generally, the main concern of these models is 



m 



and, 



whether there arises an order- disorder phase transition of the system dynamics 
if yes, how the time needed to attain consensus depends on the system size jsl, the initial 



configuration 



ion 0, y, 3 



or the topology of the underlying network 7|,l9||. 



Very recently, some extensions of the majority rule model fMRM) were proposed and 
investigated within different scenarios 10, llll, ll2|] . In Ref. 10|], Lambiotte studied a variant 
of the MRM on heterogeneous networks, namely dichotomous networks, which are composed 
of two kinds of nodes characterized by their distinct link degrees, ki and k2. It was found 
that the degree heterogeneity (characterized by the ratio 7 = ki/k2) affects the location 
of the order-disorder phase transition point and that the system exhibits non-equipartition 
of average opinion between the two kinds of nodes. The effects of community structure on 
opinion spreading and formation in the framework of MRM was also considered by Lambiotte 



and coworkers 



Motivated by the fact that a social system is inhomogeneous in many 



aspects of its inherent nature, Guan et al. introduced two ty pes of agents in the MRM where 
one type of agents have less ability to persuade the other 12|. It was shown that, as the 
inhomogeneous effect is strengthened, the location of phase transition point is shifted along 
the direction where the ordered state is more difficult to be realized. 

In the present paper, we continue the research in line of Ref. 12| to study the effects 
of free will and massive opinion of multi-agents in the MRM. The motivation comes from 
the following observations. In some cases, our opinions are strongly influenced by our social 
surroundings (the opinions of our friends, colleagues and neighbors), e.g., the fact that a 
majority of our acquaintances are smokers will easily convince ourselves to smoke; the fact 
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that a large number of our friends have MySpace will likely urge us to apply for one account. 
This phenomenon can be modelled by implementing a local majority rule when the agents 
update their opinions. In other cases, however, we make decisions in a way less dependent 
on the others, such as the type of coffee we like to drink, the fashion of clothes we like to 
wear, etc., where our taste, or "free will", dominates. In some other cases, not only the local 
majority opinion but also the global one would affect much of our behavior. Therefore, it is 
important to investigate how the free will and massive opinion of multi-agents influence the 
processes of their opinion spreading and formation. 

To proceed, we first introduce our model of opinion dynamics. The population is com- 
posed of agents located on a fully connected network. This simple setting allows us to 
look for a mean-field analytical solution for the problem. Initially, two types of opinions, 5"! 
and 5*2, are assigned to the agents with an equal probability. The densities of populations 
holding opinions Si and 5*2 are denoted by pi and p2, respectively. We investigate two types 
of models (model I and model II) involving two different opinion-updating dynamics. At 
each time step, one agent is randomly selected to update its opinion state. Following Refs. 



[lO 
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12| . two types of processes may take place: (i) With probability (1 — q), the selected 



agent interacts with two neighboring agents randomly selected from its neighborhood, and 



the three agents adopt the opinion of the local majority. The magnitude o 



acterizes the frequency of confrontation, or the strength of aggressiveness [s]. This step is 
compulsory for both models I and II. (ii) With probability q, the selected agent picks an 
opinion by using one of the following two different strategies. 

In model I, we assume that where there is no neighboring-interaction happening, an agent 
updates its opinion state according to its own will. In particular, the selected agent may 
change its current opinion to the opposite one with a free will whose strength is weighed by 
a G [0, 1]. Smaller values of a implies stronger confidence for the agent to stick to its current 
opinion. In model II, we assume that during updating, the selected agent has a tendency 
to adopt the global massive opinion with a probability proportional to its corresponding 
population density in the system: 

P{Si) - 



[1 — q) char- 



P1+P2 



where /3 characterizes the strength of the massive opinion. Larger values of f3 correspond to 
stronger tendency to become the majority. Hereafter, we always denote the global massive 



opinion by 5*1 . When there exist only neighboring interactions, i.e., in the case of g = 0, it 
is aheady known that the system asymptotically reaches global consensus where all agents 
share the same opinion j^. In the following, we investigate how the phase diagram of the 
system varies when the free will and massive opinion of agents are taken into account in the 
case of g 7^ 0. In this case, the quantity q serves as the control parameter in the dynamical 
process over the network. 

For model I, it is easy to write down the following mean-field rate equation for the system: 

At+i = At + {I - q)W + q[a{l - a) - aa], (1) 

where At is the average number of agents with opinion 5*1 at time t, a = At /N is the 
corresponding average proportion of agents with this opinion, and W is the total contribution 
to the evolution of At due to neighboring interactions. The term proportional to (1 — g) 
accounts for local majorities and the last term, for the flips with free will. The probability 
for two agents with opinions 5'i(S'2) and one with S2{Si) to be selected is 3a^(l — a) and 
3a(l — a)^, respectively, so that 

W = 3a^{l - a) - 3a(l - af = -3a(l - 3a + 2a^). (2) 

So, the evolution equation for At becomes 

At+i = ^ + (1 - g)[-3a(l - 3a + 2a^)] + qa{l - 2a). (3) 

Due to the existence of symmetry for the two opinions, it is easy to see that a = 1/2 
is always a trivial stationary solution of the above equation. In order to find nontrivial 
solutions of the system, we use the method proposed in [l^. Instead of considering directly 
the quantity At, we rewrite Eq. as follows and consider the quantities A= A — N/2 and 
5 = a- 1/2: 

Ai+i = Ai +-[3 - (3 + Aa)q - 12(1 - q)S'']. (4) 

We can easily verify that the symmetric solution a = 1/2 ceases to be stable when q < 
3/(3 + Aa), and in that case the system reaches the following asymmetric solutions: 

1 / 3-(3 + 4a)g 
^^■^ = "^'^ = 2 ^ V 12(1 -g) • 
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For convenience, we define the order parameter of tlie system by tlie difference between the 
average density of the majority and that of the minority, in the steady state, as follows: 



3-(3 + 4a)^ 
"''"^^V 3(1-,) ■ 

Thus, the system undergoes an order-disorder transition at the critical point 

Below this value, one type of opinion dominates the other, or an ordered collective phe- 
nomenon emerges due to the interactions between neighboring agents. According to Eq. 
([6]), the location of the order- disorder phase transition point depends strongly on the value 
of a, i.e., the strength of the free will of the agents to change their current opinion states. 



The stronger the inclination (to make a change) is, the more difficult an ordered col 



ective 



behavior emerges. Evidently, model I is equivalent to the MRM considered in 10| with 
a = 0.5, and in that case, Qc = 3/5. 

To test the above analysis, we perform computer simulations on model I and compare the 
simulation results with the analytical estimations in Fig. [H Simulations were carried out 
for a population of = 5 x 10'^ agents located on the sites of a fully connected network. We 




FIG. 1: The order parameter pi—p2 versus the control parameter q for several values of a. The solid 
lines are the analytical solutions obtained by calculating Eq. The symbols are the simulation 
results obtained from averaging over twenty independent experiments within a population of size 
iV = 5 X 10^. 



5 




0.0 0.2 0.4 0.6 0.8 1.0 



q 

FIG. 2: The order parameter pi —p2 versus the control parameter q for several values of /?. The solid 
lines are the analytical solutions obtained by calculating Eq. (jlip . The symbols are the simulation 
results obtained from averaging over twenty independent experiments within a population of size 
AT = 5 X 10^. 

study the order parameter pi — p2 as a function of the control parameter q. Initially, the two 
opinions 5*1 and 5*2 are randomly distributed among the agents with equal probability 1/2. 
After evolution, the system reaches a dynamic equilibrium state where the densities of the 
opinions fluctuate stably. In our simulations, one Monte Carlo step is accomplished after all 
agents have updated their opinions. The simulation results were obtained by averaging over 
the last 10^ Monte Carlo steps out of 10^. From Fig. [H we can see that the simulation results 
are in very good agreement with the analytical solutions. The small differences between the 
simulation results and the analytical predictions are due to the finite-size effect. 

Now, we consider model II, for which the mean-field rate equation for At is given by 

(1 - a)a^ - a(l - 

where the formulation of W is the same as that in model I. As assumed above, the global 
majority opinion has a greater attraction than the minority one, therefore we restrict our 
attention to the region of /3 G [0, 1) [3]. The last term in Eq. accounts for the flips 
influenced by the massive opinion, i.e., the global majority opinion of the population. Since 
the average density of the majority opinion a is less than unity (except for the full consensus 
case where a = 1), we can approximately solve Eq. (jH]) by using of the following truncated 
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Taylor expansions: 

(1 -a)^ = 1 -/5a + 0(a2), 

= 1- /3{l-a) + 0{a'^). (9) 

Substituting (j2]) and Q into (jH]), with some algebraic calculations, we get 

At+i = At + {l- g)[-3a(l - 3a + 2a')] + 97(1 - 2a), (10) 

where 7 = (1 - /5)/(2 - /?). 

The remaining part is the same as what has been done for model I. We hence omit the 
lengthy calculations and simply write out the final result for the order parameter of model 
II: 



Pl - P2 



/3 -(3 + 47)g 
3(1 -g) 



3-[3 + 4(l-/5)/(2-^^ ^^^^ 



3(1 -g) 

Thus, in model II, the system undergoes an order-disorder transition at the critical point 

3(2-/?) +4(1-/?) 10-7/?- ^ ^ 



0.9 



O (if 
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□ N = 5000 
• N = 10000 
A N = 20000 
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q 



FIG. 3: The order parameter pi — p2 versus the control parameter q for different system sizes A'^. 
The parameter /3 = 0.9 is fixed. The arrows point to the locations of phase transition obtained by 
simulations. 
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For f3 = 0, our model reduces to the MRM studied in lOj, so it is not surprising to recover 
the known result of gc = 3/5. While for other positive values of /5 less than unity, we obtain 
a threshold of q (at which the order-disorder transition takes place) greater than 3/5. That 
is, if the information of global massive opinion is accessible to the agents who have tendency 
to becoming the majority, it is easier for the whole system to attain an ordered state. 

The simulation results of the order parameter pi—p2 as a function of the control parameter 
q for several values of (3 are shown in Fig. [2l The analytical solutions of Eq. f|TT]) are also 
plotted there for comparison. Once again, they well match with each other. The differences 
between the simulation and analytical results come from two aspects: one is the finite-size 
effect, and the other is the approximation of Eq. From Fig. HI we also notice that for 
any large enough value of /3, there arises another phase transition where a full consensus 
state emerges as the parameter q goes below some critical value g^, i.e., one of the two 
opinions dominates the whole system for any q < q'^. 

We argue that the fact that the considered system is of finite size contributes to the phase 
transition. In fact, reviewing Eq. ([H]), we find that At+i will be going to iff g satisfies the 
following condition: 

3a(l - a)(2a - 1) 



^ < = 77Z (13) 

a/3+{l-a)/3 



3a(l - a){2a - 1) ^ {l-a)aP-ail-a) 



Note that a is the average density of the majority, a G [0.5,1]. In this region of a, it is 
easy to verify that /(a) is a monotonically decreasing function for positive (3. Thus, as long 
as q is smaller than the minimal value fmin of /(a), the whole system will evolve to a full 
consensus state. 

From an ecological point of view, the number of species of a population must be kept above 
a minimum level to prevent going to extinction due to random fluctuations of environmental 
conditions. Inspired by this viewpoint, for the present MRM, we judge that the number of 
the minority should not be smaller than 2, below which they are doomed to extinct (since two 
agents of the same opinion is guaranteed to convince a neighbor), but above which they have 
a chance (though very small) to survive. Thus, a reasonable estimation is fmin = /(1~2/A^). 
In Fig. [31 we present simulation results for three systems with different sizes, = 5 x 10^, 
10^, and 2 x 10^, respectively. We can see that the critical value of q'^ is decreasing with the 
increase of the system size. For P = 0.9, the simulations yielded that q'^ = 0.73, 0.70, and 
0.66 for = 5 X 10'^, 10^, and 2 x 10^, respectively, which are in good accordance with the 



estimations of /(I — 2/N) = 0.717, 0.692, and 0.665 for the three corresponding values of 
N. 

To sum up, we have investigated the effects of free will and massive opinion of multi- 
agents in a network based on the majority rule model, which is a simple yet useful statistical 
model for studying the emergence of collective behaviors. Two types of models have been 
considered. In model I, the agents have a free will to either adopt the opposite opinion or 
remain its current opinion, whereas in model II, besides the local majority, the agents also 
have a tendency to becoming the global majority when there are no neighboring interactions 
among them. We found that the location of the order-disorder phase transition point strongly 
depends on the involved dynamics, which may give rise to either smooth or harsh conditions 
to achieve an ordered state. In addition, for finite system sizes, we found that there exists 
a threshold below which a full consensus can be attained if the agents have a inclination 
to becoming the global majority. Our analytical estimations are in good agreement with 
simulation results. In the present work, however, we only studied the majority rule model 
over fully connected networks. To model the real world in a more accurate fashion, future 
research may take into account the coevolution of the majority process and the underlying 
interaction pattern, i.e., from the viewpoint of adaptive and co-evolutionary networks 

This work was supported by the NSFC-HK Joint Research Scheme under the grant N- 
CityU 107/07. 
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